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This paper, based on a presentation at the Spintronics 2001 conference, provides a review of our 
studies on II- VI and III-V Mn-doped Diluted Magnetic Semiconductors. We use simple models 
appropriate for the low carrier density (insulating) regime, although we believe that some of the 
unusual features of the magnetization curves should qualitatively be present at larger dopings (metal- 
lic regime) as well. Positional disorder of the magnetic impurities inside the host semiconductor is 
shown to have observable consequences for the shape of the magnetization curve. Below the critical 
temperature the magnetization is spatially inhomogeneous, leading to very unusual temperature 
dependence of the average magnetization as well as specific heat. Disorder is also found to enhance 
the ferromagnetic transition temperature. Unusual spin and charge transport is implied. 



I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) are com- 
prised of an inert host semiconductor doped with both lo- 
calized spins and carriers (electrons or holes) that are ei- 
ther itinerant, or localized on a much longer length scale. 
In that sense, they belong to the general family of cor- 
related electron systems, which include a number of fas- 
cinating materials such as cuprates, manganites, heavy 
fermions and other Kondo lattice systems. 

Electronic materials containing local moments have 
been studied for some time. What makes the DMS so 
fascinating is that they belong to a regime that has previ- 
ously been neglected. While the name diluted magnetic 
semiconductors implies (correctly) that the system has 
only a small percentage of localized spins, they are at 
the opposite extreme of the dilute magnetic alloys such 
as Fe or Mn in Cu, the canonical systems involving itiner- 
ant fermion and localized spin degrees of freedom, which 
have been studied extensively.EJ In the dilute magnetic 
metallic alloys, the low density of spins are a perturba- 
tion on the Fermi liquid representing the non-magnetic 
host metal, so depending on the concentration of the local 
moments, they may be studied in terms of dilute Kondo 
systems, or amorphous magnetic systems with a spin-spin 
coupling mediated by the Fermi sea of conduction elec- 
trons (RKKY coupling), which lead often to spin glass 
behavior.a 

By contrast, in the regime of interest, the carrier den- 
sity in DMS is significantly lower than the (low) localized 
moment density, so the spins become an integral part of 
the description of the system and its magnetic phase, 
rather than a mere perturbation on a metallic Fermi 
sea. In that sense, the situation is even more extreme 
than e.g., in Kondo lattice and heavy Fermion materials, 
where the two species have comparable densities. This 
large, inverted, ratio of local moments to carriers is in 
fact similar to that in the high Tc cuprates. However, 
unlike the cuprates, the density of local moments is low 



and incommensurate with the lattice, and the carriers 
and the spins are not in the same band. As a conse- 
quence of the low moment density, the exchange between 
local moments is not standard direct or superexchange, 
as in the cuprates, but is mediated by the carriers, even 
though their density is so small. Thus, the DMS are in 
rather different region of phase space of electronic mate- 
rials with local moments, than other correlated electron 
systems. 

Despite this difference, most models of diluted mag- 
netic semiconductors start from the high carrier density 
limit, where the carriers may be modeled as-|free carriers 
moving in the conduction or valence band.tl This is un- 
derstandable, since in the high density limit the carrier 
kinetic energy is the largest energy in the problem, and 
calculations may be done perturbatively starting from 
the non-interacting Fermi gas. However, most of the in- 
teresting behavior is seen at low carrier densities, where 
the system is insulating, or not too far from the metal- 
insulator transition. Consequently, we have concentrated 
in this work on the low density regime, starting from 
bound carriers, and moving on to carriers in an impurity 
band formed from the bound impurity states. 

As we wish to cover the case of insulating behavior 
at arbitrary filling factor i.e., away from the half filled 
impurity band case (one carrier per site), disorder has 
to be included at the outset. In particular, we model 
the system with randomly distributed dopants, as in 
the experimental system, since it has been recognized 
that the random distribution is essential to understand 
the magnetic properties^of conventional, non-magnetic 
doped semiconductors.Qi3 Such models exhibit both in- 
sulating and metallic phases, and with the random dis- 
tribution of impurity sites included, are in a position to 
reveal the effect of disorder in low carrier density sys- 
tems. In the case of the predominantly antiferromag- 
netic couplings between hydrogenic centers in conven- 
tional doped semiconductors, the randomness is found 
to suppress magnetic order below measurable temper- 



1 



atures ( ~ millikelvin ), and possibly to zero. In the 
case of ftMS, where interactions lead to ferromagne 
ordering,Bi3 in agreement with experimental findings| 
we find that randomness leads to unusual behavior in 
the magnetic response, and effects of randomness are ex- 
pected in the transport behavior as well. 

In this paper, we review the results of our approach 
to DMS based on both II- VI semiconductors (like CdTe 
or ZnSe), and on III-V semiconductors (such as GaAs 
or GaP), and compare the two families of DMS systems. 
The paper is organized as follows. In Section || we briefly 
review the properties of conventional (insulating) ferro- 
magnets. The results presented serve as a reference with 
which to contrast the results obtained in the rest of the 



paper. Section |II| addresses the II- VI based DMS, while 
Section [V deals with III-V DMS. In each case, we intro- 



duce the Hamiltonians we use to model these systems. 
Results obtained within mean-field approximation and 
with Monte Carlo simulations are presented. The effect 
of positional disorder of the Mn dopants is studied, as are 
the similarities and differences between II- VI and III-V 
DMS. Finally, Section ^ summarizes our results and con- 
clusions. It also includes a discussion of important issues 
such as robustness of models, relevance of disorder on 
spin scattering, and key experiments which could help 
provide a better understanding of these materials. 



II. CONVENTIONAL FERROMAGNETIC 
SYSTEMS 

A typical model of a uniform ferromagnet consists of 
a collection of identical magnetic spins of magnitude S, 
placed on an ordered Bravais lattice. While the generic 
case may be anisotropic in spin space due to a variety 
of reasons (spin-orbit couphng, crystal fields etc.) we 
consider here the simplest isotropic case where the spin 
interactions are well described by a Heisenberg Hamilto- 
nian 
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Si ■ Sj — H 



(1) 



Due to translational invariance, the exchange integral 
depends only on the distance Ri — Rj between spins. In 
insulating materials, this dependence is due to overlap 
between the electronic orbitals involved in creating the 
spin S (through Hund's rule), leading to an exponential 
decay of with increasing distance. In typical magnetic 
atoms, these orbitals are of d or / type, and they are 
localized within ^1 — 2 Aof the nucleus. As a result, 
it is customary to restrict the first sum in Eq. (|l|) to 
only nearest-neighbor spins. The external magnetic field 
H = HBz breaks the rotational symmetry, leading to the 
appearance of a non- vanishing expectation value (Sf) at 
each site. Translational invariance implies that (S^) = 
{S) is independent of the position Ri of the spin. 



While an exact solution for the Heisenberg Hamilto- 
nian (Q) is known only in one dimension, it has been 
found that the Weiss (mean-field) approximation pro- 
vides a qualitatively good understanding of the proper- 
ties of these systems. The mean-field factorization Si ■ 
S, ^ S,-{S,) + {S,)-S,^{S,)-{S,) = (S) {Sf + S^)-{S)^ 
allows for a solution of the problem in terms of an ef- 
fective magnetic field H{i) ^ H + J{S) /{gfis), where 
J = '^j^i Jij- [If only nearest- neighbor interactions are 
kept, J = zJqi, where z is the coordination number 
of the Bravais lattice and Joi is the exchange integral 
for nearest-neighbor spins]. In the absence of an exter- 
nal magnetic field, a non-vanishing solution for {S) is 
found for T < Tc ^ JS{S + l)/3fcB. In other words, 
the system is ferromagnetically aligned below the criti- 
cal temperature Tq, and the spontaneous magnetization 
(M) — g^B{Si) increases rapidly (Fig. |l|) with decreas- 
ing temperature and is already close to the saturation 
value Mo — g^sS below T < 0.5Tc- Concurrently, the 
specific heat has a peaked structure around Tc and drops 
rapidly to zero for T < 0.5Tc reflecting the fact that the 
only accessible degrees of freedom for low T are the long- 
wavelength (collective) spin-wave exoitations which have 
a restricted phase space (see Fig. 

It is well-known that mean-field approximations over- 
estimate the strength of the correlations, leading to 
rather high estimates for the Curie temperatures Tc- De- 
tailed studies of these Hamiltonians with Monte-Carlo 
simulations, which properly account for the effects of 
thermal fluctuations, find quantitative changes of up to 
a factor of 2 in the value of Tc- However, as suggested 
in Fig. 1^, the qualitative features of the magnetization, 
specific heat and susceptibility curves remain as in the 
Weiss mean-field treatment, in good agreement with ex- 
perimental measurements. 



III. FERROMAGNETIC SYSTEMS: II- VI DMS 

A. The model 

The II- VI DMS are based on semiconductors AB, 
where A is a group-II element and B is a group- VI ele- 
ment (such as CdTe or ZnSe). In the II- VI DMS, some of 
the divalent sites (Cd/Zn) are substituted by a magnetic 
element, typically Mn. This fraction is denoted by x, so 
the DMS we consider is Ai_a;Mna;B. Mn is also a group- 
II element, but in addition it has a half-filled 3d shell, 
with a total spin given by Hund's rule: S = 5/2. In the 
absence of other types of dopants, the system Ai-a^Mn^^B 
is an insulator which exhibits antiferromagnetic tenden- 
cies at low temperatures. This is seen, for instance, from 
measurements of the susceptibility which is found to de- 
pend on temperature as x{T},(^l/{T + Tn), with a Neel 
temperature of a few kelvin.cl'llj The origin of this anti- 
ferromagnetic (AFM) tendency is the (expected) antifer- 
romagnetic exchange between the Mn spins. However, 
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for low doping concentrations x, the average distance be- 
tween Mn spins is large and this AFM direct exchange is 
rather small. 

When a low density of charged dopants, such as group- 
V Phosphorus (P) substituting for the group- VI element, 
is introduced in the system, each of them binds a hole 
(or electron) in a shallow hydrogenic Is state 0(r) ~ 
exp(— r/as), characterized by a Bohr radius as ~ 10—20 
A. Exchange interactions arise between the spins of these 
charge carriers ajid the Mn spins, and are described by 
the HamiltonianQ: 

7J = ^J(rl,i?,)s,.5,. (2) 

Here, Sj is the spin of the Mn at position Rj and Si is 
the spin of the electron/hole centered at n. The exchange 
interaction J(ri, Rj) is dependent on the overlap between 
the orbital <j){f—ri) of the charge carrier and the orbitals 
4>d{r—Rj) of the 3d electrons responsible for the Mn spin. 
Since these 3d orbitals are localized on a scale of a few 
A around the Mn nucleus, the exchange is proportional 
to the carrier charge density at the Mn site, i.e. 

J(r„Rj) = Jo|0(i?j - rl)|2 = Joe-2|r.-R.I/-B^ (3) 

where Jq characterizes the strength of the exchange. 
Typically, for electrons Jq < 0, while for holes Jq > 0. 
However, since in the following we treat the spins as clas- 
sical variables, the sign is irrelevant. For specificity, in 
the rest of the paper we assume Jq > corresponding to 
holes as charge carriers. 

The Hamiltonian neglects the direct AFM inter- 
actions between the Mn spins. For low values of the 
fraction x, it can be simply accounted for in the follow- 
ing manner. For Mn spins which are very close to one 
another (such as nearest neighbors), the direct AFM ex- 
change is the dominant (large) interaction, and leads to 
the formation of a singlet state. This singlet becomes in- 
ert as far as magnetic interactions are concerned. For Mn 
spins which are fairly far apart from other Mn spins, the 
dominant magnetic interaction is the exchange with the 
charge carrier spins. As a result, to first order the Hamil- 
tonian (^ accounts for both types of interactions if we 
restrict the summation over the Mn spins to only those 
Mn spins which are not part of a spin-singlet. At low x, 
this includes a large majority of Mn spins. If the fraction 
X of Mn becomes too large, both types of interactions will 
be of comparable size for all the Mn spins, and therefore 
this separation is no longer possible. In this case, the 
frustration imposed by the competing exchanges leads to 
the appearance of a spin-glass state, which has been ob- 
served experimentally for x > 0.2.t3 In the following, we 
restrict ourselves to the low x {x < 0.1) limit. 

Simple thermodynamic considerations show that, qual- 
itatively, at a temperature fc^T < J{r) [see Eq. (||)], 
all Mn spins within distance rx ~ (03/2) In ( Jq/ZcbT) 
of a dopant order their spins antiferromagnetically with 



respect to the dopant hole spin. As a result, a region 
with a large magnetization (from all the parallel polar- 
ized Mn spins) appears near the dopant. This is known as 
a Bound Magnetic Polaron (BMP ),□ whose radius rx (see 
above) increases logarithmically with decreasing temper- 
ature. As a result, one expects that long-range ferromag- 
netic order appears in the system for temperatures low 
enough that a continuous percolating network of BMPs is 
formed (as shown schematically in Fig. provided that 
nearby BMPs prefer to orient ferromagnetically with re- 
spect to one another. At first sight, this seems to not be 
the case, since direct exchange between the charge carri- 
ers localizaiitt hydrogenic orbitals has an antiferromag- 
netic sign.QO'Ed However, this antiferromagnetic coupling 
is overwhelmed by effectively ferromagnetic interactions 
between BMPs coming in part from Mn spins in between 
the piolarons which favor ferromagnetic alignment of Mn 
spinscl and partly from the modification of the effective 
direct exchange as a result of the local field due to the 
polarized Mn.l3 

These mechanisms favoring parallel orientation of the 
BMPs at low temperatures are rather weak, and as a re- 
sult the Curie temperature below which long-range fer- 
romagnetism is observed in these systems is very low, 
to our knowledge below 5 K for all II- VI DMS studied 
so far. Moreover, as the transition is of a percolation 
type, and the percolation fraction is ^ 20% for three di- 
mensions, this implies that just below Tc about 80% of 
the Mn spins do not participate in the ferromagnetism. 
These are the Mn spins which are outside the percolated 
cluster, i.e. far from the charged dopants (see Fig. ||). 
They are very weakly interacting (essentially disordered) 
unless the temperature becomes so low that a nearby 
BMPs grows large enough to include them. This results 
in a very unusual FM phase, in which a substantial part 
of the spin entropy survives down to very low T. 

B. Monte Carlo simulations 

We performed Monte Carlo simulations on the Hamil- 
tonian (^) , to study this unusual FM phasp-treating both 
Mn and carrier spins as classical variables .EJ This appears 
to be a reasonable approximation, since = 5/2 is a large 
spin and the Mn spins dominate the magnetic response. 
Simulations were carried out for zinc-blende lattices with 
lattice constant a = 5 A, for Mn concentration x = 0.001, 
dopant density Ud = 10^* cm~^ and as = 20 A. The ex- 
change Jo defines the unit of energy. With these parame- 
ters, the Mn concentration umu — ^x/a? is 32 times the 
dopant concentration. Nevertheless, the magnetic cou- 
pling is mediated by the latter because of the large Bohr 
radius, as required for the polaron picture to hold. 

The magnetization curves obtained have unusual, con- 
cave upward shapes (see Fig. ^ left panel), very unlike 
the typical magnetization curve of Fig. |l|. For these pa- 
rameters, the critical temperature Tc = 0.014 Jq is found 
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using finite size scalingB We find that the magnetization 
reaches its saturation value only at exponentially small 
temperatures, reflecting the existence of the quasi-free 
Mn spins outside the percolated (magnetically ordered) 
region. 

The specific heat of the classical Heisenberg model has 
the unphysical limit Cy NkB as T — > (empty squares 
in Fig. 1^). While this agrees with the equipartition 
theorem, it implies that quantum mechanics (with dis- 
crete energy levels) is needed to capture the correct limit 
Cy — > as T 0. One way to mimic the discretiza- 
tion, but avoid the complexities of the quantum Monte 
Carlo treatment, is to use a discrete (classical) vector 
model, in which each Mn spin can only be oriented along 
one of the six [100] directions. An efficient Monte Carlo 
method for this discrete model is described in Ref. 
While the magnetization curves are very similar to the 
ones obtained in the continuous spin Heisenberg model 
(see Fig. 0), the specific heat results are very different 
(see Fig. As expected, for the discrete model Cy 
as T ^ 0. However, unlike in the case of a typical FM, 
the peak in Cy is not near Tc, but at temperatures well 
below Tc- This refiects the residual entropy of the free 
Mn spins outside the percolated region. 



C. Effect of disorder 

In H-VI DMS there are two sources of positional dis- 
order: disorder in the positions of the Mn spins and dis- 
order in the position of the charged dopants. In the limit 
when there are many Mn ions per dopant, the Mn spin 
disorder is not expected to have a significant effect on the 
magnetization curves, or the critical temperature. The 
reason is that at the very low temperatures where per- 
colation appears, the radius of each BMP is significantly 
larger than as, favoring interactions with a large num- 
ber of Mn spins. Disorder in the Mn positions will lead 
to some fluctuations in the average number of Mn spins 
found in each BMP, but this should have a relatively 
small effect. 

On the other hand, disorder in the position of the 
charged carriers (centers of the BMPs) has a large ef- 
fect on the critical temperatures. As seen from Fig. H, 
disorder in the positions of the BMPs facilitates the ap- 
pearance of a large percolated cluster for smaller BMPs 
sizes (larger temperatures), since only a subset of the 
BMPs must percolate in order for ferromagnetic order to 
appear in the system. On the other hand, the ordered 
BMP lattice only percolates when each and every BMP is 
included. This obviously happens when a larger fraction 
of the space is filled by BMPs, i.e. at a lower temperature. 
However, we emphasize again that even for the ordered 
BMP lattice, a significant volume containing a large frac- 
tion of the Mn spins is still outside the percolated volume 
(in the interstitial spaces) and therefore the phenomenol- 
ogy related to the existence of weakly-interacting spins 



down to exponentially low temperatures is still valid. 

We have verified, using Monte Carlo simulations, that 
the critical temperature of a system in which the charge 
carriers are placed in an ordered supcrlattice is lower than 
that of a sample with random positions of the charge car- 
riers, when all other parameters are identical. For the 
case investigaied, the relative increase of Tc with disor- 
der was 50%Ji3 However, this relative increase is expected 
to depend on the various parameters of the problem. 



IV. FERROMAGNETIC SYSTEMS COUPLED TO 
FERMIONS: III-V DMS 

A. Introduction 

When Mn is doped in a Ill-V semiconductor, such as 
GaAs, the major difference with respect to the II- VI 
DMS is that the Mn atom provides both the S = 5/2 
spin and the dopant charge carrier (a hole, since diva- 
lent Mn substitutes for trivalent Ga). While this implies 
nominally equal concentrations of holes and Mn spins, 
experimentally it is found that the hole concentLation 
is only p = 10 - 30% of the Mn conccntration.&y The 
compensation process (es) responsible for the removal of 
such a large fraction (~ 70 — 90%) of the holes from the 
carrier band are not fully understood, but it is believed 
that an important role is played by As antisite defects. 
Such defects are created when group-V As substitutes 
for group-III Ga, and removes two holes introduced by 
Mn impurities, thus effectively decreasing the hole con- 
centration. Compensation is responsible not only for the 
substantial decrease of the hole concentration, but also 
leads to the appearance of charged compensation cen- 
ters (e.g. As^+ for As antisites). The Coulomb potential 
created by these charged compensation centers may also 
play a role in the physics of these systems, as we discuss 
in the following. 

As in the (II,Mn)VI systems, the main magnetic in- 
teraction in the (III,Mn)V DMS is the exchange be- 
tween the Mn spins ancLthe hole spins, which is known 
to be antiferromagneticJa Assuming, again, very sharply 
peaked Mn 3d orbitals, this exchange is proportional 
to the probability of finding the charge carrier at the 
Mn site. This probability is extracted from the wave- 
functions of the orbitals occupied by the hole charge car- 
riers. The appropriate framework to describe the hole 
states depends on their concentration. At low hole con- 
centrations, screening processes are ineffective. The un- 
screened Coulomb potentials of the Mn dopants are re- 
sponsible for the splitting of hydrogen- like impurity levels 
from the top of the valence band, and the holes occupy 
these impurity states. In the limit of high hole concen- 
trations when the carrier kinetic energy is the largest 
energy in the problem, the Coulomb potential of the Mn 
dopants effectively gets reduced because of screening. As 
a result, the holes occupy a Fermi sea at the top of the 
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valence band. Qualitatively, it is apparent that the two 
situations could lead to quite different physics. Holes oc- 
cupying Bloch states in the valence band are found with 
equal probability anywhere inside the host semiconduc- 
tor, and therefore one expects the system to be rather 
homogeneous. On the other hand, holes occupying im- 
purity states are found with high probability near the 
Mn sites. As a result, we expect a rather inhomogeneous 
distribution of the holes in the host semiconductor, and 
the positional disorder of the Mn dopants may play an 
important role, since it defines the length-scale for these 
inhomogencities . 

Gai_j;Mn2,As has a Metal-Insulator Transition for 
X ~ 0.03, and shows re-entrant insulating behavior for 
X > 0.07.Q In the insulating regimes, the low-temperature 
conductivity. is consistent with Mott long-range variable 
hopping ,E£lEZl suggesting the existence of impurity-like 
levels. Even for the most metallic sample {x = 0.053) 
the screening length (/ ~ 10 A) as evaluated from the 
Thomas-Fermi theory is of comparable size, not much 
smaller jihan the Bohr radius of the impurity level {as ^ 

8 A)Mki 

B. The model 

Motivated by these observations, we have attempted 
to understand the low x regime within a model based on 
the existence of impurity hydrogen-like orbitals at each 
Mn site. While this is similar to our approach to the II- 
VI DMS systems, one difference is that since the number 
of holes is smaller than the number of Mn, there must 
be a mechanism to allow the holes to "choose" the Mn 
dopants near which to stay. Such a mechanism is nat- 
urally provided by hopping processes facilitated by the 
overlap between impurity wave- functions centered at dif- 
ferent Mn sites. Therefore, the Hamiltonian describing 
such a system is of the form 

- .g/isif ^ |4c,. - ~g^iBHj2S'{^)- (4) 

i i 

Here, i indexes different Mn positions Ri, and cj^ is the 
creation operator for a hole with spin a in the impurity 
level centered at Ri, while Si is the spin of the corre- 
sponding Mn dopant. 

The first line in Eq. (|^) is the Hamiltonian of the 
charge carriers. The first term describes hopping of 
holes between impurity levels. For simplicity, we assume 
again Is impurity states with 0(r) — e:Kp{—r/aB)- In 



fact, the hole impurity wave- function is more compli- 
cated, due to the band-structure of the valence band 
from which it splits (for details, see Ref. |l|). For 
the hopping integral we use the simple parameterization 
tij ~ 2(1 + r/ae) exp (— r/ae) Ry, where r — |Ri — Rj|, 
appropriate for hopping between two isolatedJs impu- 
rities which are not too close to one another.EEl For Mn 
doped into GaAs, the Bohr radius is as = 7.8 A and the 
bindipSijenergy which defines the Rydberg is 1 Ry =110 
meV.llflo We have investigated other |parameterizations 
for the hopping matrix t{r) elsewherejiil and found that 
while they lead to quantitative changes, qualitatively the 
results are similar. 

The second term describes an on-site potential u{i) due 
to the Coulomb potential of the other Mn impurities, as 
well as other nearby charged compensation centers. An 
on-site Coulomb repulsion U of the Hubbard type may 
be added to describe the electron-electron repulsion be- 
tween electrons occupying the same impurity orbitals. 
For isolated Is impurities, t/ « 1 Ry. However, depend- 
ing on the effectiveness of screening, the electron-electron 
interactions may be longer-range. A fully self-consistent 
treatment of this problem should involve a proper de- 
scription of the screening processes, and would allow a de- 
tailed computation of the strength of the hopping matrix, 
the on-site Coulomb potential and the electron-electron 
interactions. However, since the full self-consistent de- 
scription is extremely difficult to achieve, especially as 
details about compensation processes are still not clari- 
fied, we use the simplified assumptions described above. 
We believe that they should provide a good qualitative 
description of the properties of these compounds, and 
with proper fitting of various energy and length scales 
may even lead to a quantitative description. 

The second line of Hamiltonian (^ describes the AFM 
exchange between the Mn spin Sj and the hole spin 

Si = ^cl^<JapCii3 [(? are the Pauli spin matrices]. As in 
II- VI DMS, the AFM exchange is proportional to the 
probability of finding the hole trapped at Ri near the 
Mn spin at Rj, so = J exp (^—2\Ri — Rj\/aB^- Based 

on calculationsEl of the isolated Mn impurity in GaAs, 
we estimate the exchange coupling between a hole and 
the trapping Mn {Ri — Rj) to be J = 15 meV. 

The second term describes the direct Mn-Mn exchange 
in the semiconductor host, which is expected to be short 
range, and consequently not important at low x when Mn 
are a few sites away from each other. We have therefore 
omitted this term (i.e. set Kij = 0); however, for higher 
concentrations this may be important. Finally, the third 
line in Hamiltonian (^) describes the interaction with an 
external magnetic field. 

Given the large number of terms in the Hamiltonian, it 
is useful to try to understand the effect of each. To begin 
with, we neglect the random on-site potential {u{i) = 0), 
the electron-electron interaction (?7 = 0), the direct Mn- 
Mn AFM interactions {Kij = 0) and turn off the external 
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magnetic field {H = 0) (we will discuss the effects of these 
various terms later on). As a result, the Hamiltonian con- 
tains only its two main terms {tij and Jy), describing the 
dynamics of the charge carriers and the AFM interaction 
between the Mn spins and the charge carrier spins. 

C. Similarities and differences between II- VI and 
III-V DMS 

We investigated the Hamiltonian— using both the 
mean-field approximation (MFA)EjM and Monte Carlo 
(MC) simulations.c3 Typical magnetization curves ob- 
tained with MC methods for a Mn concentration x = 0.01 
and hole concentrations p = 10 and 30% are shown in Fig. 
^. The corresponding curves obtained using the mean- 
field approximation for the similar parameters are shown 
in Fig. |j. While there are substantial quantitative differ- 
ences between the two, these are easily understandable. 
The long tail of the MC curves at high T are due to fi- 
nite sizes of the samples studied; these disappear as the 
sample size is increased. On the other hand, the criti- 
cal temperature (Tc) predicted by MFT are significantly 
higher than those obtained by MC simulations (as would 
be expected). Part of the difference in Tc between the 
two methods is actually due to the fact that the Mn spins 
in the MC simulations are taken to be classical variables, 
and quantum operators in the MFA. If we use classical 
Mn spins in MFA, we find Tc reduced by a factor of 
~ 2. The remaining reduction is presumably due to the 
usual neglect of fluctuations in MFA, which is properly 
captured in MC simulations. 

The striking feature, common to both results, is that 
the magnetization curves have unusual shapes - linear or 
concave upward. This is qualitatively similar to those 
found for the Il-VIrfl^^ (Fig. and what has been 
seen in experiments|l3'E2l but very different from the con- 
vex upward M{T) of conventional ferromagnets (Fig. |^). 
Again, as in the II- VI, the magnetization does not reach 
its saturation value until very low temperatures. Con- 
currently, the specific heat curves also exhibit a peak at 
temperatures much lower than the critical temperature, 
reflecting the entropy of the disordered spins, present in 
the system down to these low temperaturesO 

By looking at the magnetization profile around the 
Tc, long-range ferromagnetism in the disordered sam- 
ple of III-V DMS is seen to appear when a percolated 
cluster of polarized Mn spins is formed. However, un- 
like in the spin-only model based on isolated hydrogenic 
centers used for II- VI DMS in the previous section, the 
holes are delocalized within this cluster for the parame- 
ters appropriate for the III-V based DMS. Since the holes 
can more effectively minimize their kinetic energy when 
maintaining the direction of their spin during hopping, 
this delocalization of the holes within the percolated clus- 
ter provides a very effective mechanism for alignment of 
all Mn spins within the cluster in the same direction. 



This kinetic-induced alignment mechanism is much more 
effective than mechanisms of alignment of nearby BMPs 
in insulating II- VI DMS, suggesting higher critical tem- 
peratures in this case. Other reasons for enhancement 
of critical temperatures in III-V DMS include the peak- 
ing of the impurity wavefunctions at the Mn sites in this 
case where Mn is also the dopant, and the ability of car- 
riers in the compensated case to choose states with wave- 
functions peaked in the regions with higher-than-average 
Mn concentrations - the higher probability of finding the 
holes in these regions leads to enhanced effective inter- 
actions with the Mn spins. When all these factors are 
included, we find indeed that the striking differences in 
critical temperatures, by two orders of magnitude, in the 
two systems can be comfortably explained, at least within 
MFA. 



D. Effect of disorder 

Within the mean-field approximation, positional disor- 
der in the Mn spins for III-V DMS leads to a significant 
increase of the critical temperature. Mn are the charged 
dopants in this case, so the situation again appears to 
be similar to that in II- VI DMS. Typical magnetization 
curves obtained using MFA are shown in Fig. ||, for a 
doping X — 0.0093 and p = 10%. In order of increasing 
Tc, the four curves correspond to increasing disorder in 
the positions of the Mn impurities. We start with a fully 
ordered, simple cubic superlatticc of Mn impurities inside 
the host semiconductor (for this concentration, the super- 
lattice lattice constant is equal to three lattice constants 
of the underlying Ga FCC sublattice). The correspond- 
ing average spins of the Mn and charge carriers are shown 
by dashed lines in Fig. ^j. Then, we introduce positional 
disorder of the Mn ions on the underlying Ga FCC sub- 
lattice in varying amounts - (i) low-disorder where Mn 
spins are randomly place on any of the nearest-neighbor 
sites of the original superlattice sites; (n) moderate dis- 
order - where the Mn spins are allowed to occupy any 
sites on the Ga sublattice, as long as the distance be- 
tween any two Mn is larger than two lattice constants; 
and (iii) completely random positions of the Mn spins on 
the Ga FCC sublattice of the host semiconductor. 

In a fully ordered III-V DMS, below Tc each Mn spin 
is equally polarized, since translational invariance implies 
that the holes are equally distributed among the various 
Mn sites and therefore create the same effective magnetic 
field for each Mn spin. [ In this respect the ordered lat- 
tice for III-V is different from the situation encountered 
for the ordered superlattice of charged dopants in the II- 
VI DMS. In the II- VI, below Tc the Mn spins inside the 
BMPs are strongly polarized, while the Mn spins out- 
side the BMPs are practically unpolarizcd ] . The reason 
why Tc is larger in a disordered III-V DMS than an or- 
dered one, is that the hole wave-functions are pulled-in 
the regions with higher-than-average Mn concentrations. 
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where they can more effectively minimize their total en- 
ergy. The increased probability of finding the holes in this 
smaller volume occupied by the cluster leads to effectively 
larger couplings Je/ / of the Mn spins in the cluster and 
therefore increased critical temperatures. In other words, 
in the disordered III-V DMS the holes only need to po- 
larize a smaller fraction of the Mn spins in the system 
and get polarized in their own term. In an ordered Mn 
sample, the holes polarize equally all the Mn spins in the 
system, and this can only happen at rather low tempera- 
tures, given the small number of holes as compared with 
the number of Mn spins. While MFA shows a strong de- 
pendence of Tc on disorder, this is likely to be modified 
once fluctuation effects left out in MFA are included, as 
in a MC simulation. 

The unusual shape of the magnetization curves is a 
consequence of the relatively small number of charge car- 
riers as compared to the number of Mn spins. In a dis- 
ordered system, we can identify two types of Mn spins: 
strongly-interacting Mn spins from the percolated clus- 
ter, which polarize at high temperatures and lead to the 
ferromagnetic transition at Tc, and weakly-interacting 
Mn spins from the regions outside the percolated cluster. 
Since these outside regions have low hole density in our 
model, the effective coupling of their Mn spins (which is 
proportional to the probability of finding holes nearby) 
is rather small, Consequently, these spins do not polar- 
ize unless the temperature is comparable in size to their 
effective coupling. We have used this picture to obtain a 
simplified but fairly accurate description of the magnetic 
and thermodynamic pmperties of the DMS based on a 
two-component model.Lj We start from a histogram of 
the effective couplings Je// of all the Mn spins, obtained 
by averaging over many realizations of disorder. Such a 
histogram of Jeff /J obtained using Monte Carlo simu- 
lations for X = 0.01 and p — 10% is shown in Fig. ||. As 
can be seen, it is a very wide distribution, from very large 
Jeff ~ J for strongly interacting Mn spins, to extremely 
small Jeff /J ^ 10'^ values for weakly interacting Mn 
spins. Histograms obtained within the mean-field ap- 
proximation have vepjj" similar shapes, except that their 
width is even larger .Eil 

For such wide distributions, at any given tempera- 
ture ksT, we divide the spins into weakly/strongly- 
interacting categories, depending of whether their effec- 
tive coupling Jeff is smaller/larger than jksT. Then, we 
replace the complex distribution of couplings shown in 
Fig. H by two ^-functions representing the two spin com- 
ponents. The values of the nominal couplings Ji and J2 
of the weakly/strongiy-interacting spin components are 
simply the average of all the couplings of weakly/strongly 
interacting spins. The constant 7 is found from a fit of, 
for instance, the magnetization curve provided by this 
simplified model. Other thermodynamic quantities, such 
as susceptibility and specific heat are then_.shown to be 
quite well described by this simple model.o In contrast, 
we have verified that replacing all the couplings by a sin- 
gle coupling corresponding to the average over the entire 



distribution leads to curves very different than the ones 
obtained with the original distribution. We believe this 
simplified model could provide a simple tool for inter- 
pretation of experimental curves. So far, most attempts 
have been to try to fit the magnetization curves (for in- 
stance) with only one coupling. While this may recapture 
part of the curve near and below Tc!u it turns out that it 
only accounts for a rather small percentage of the total 
number of Mn spins expected to be in the system. This 
suggests that a second component is missing. In fact, fits 
in terms of two components, one ferromagnetic and one 
paramagnetic, have already been performed in ordep-jto 
explain the shapes of the measured M{H,T) curves.Ej 

In a conventional ferromagnet, an external magnetic 
field will lead to a fast increase of the magnetization 
from its value in the absence of the field, to the satura- 
tion value Mq — NgfisS, where N is the concentration 
of spins S in the system. A hysteresis curve associated 
with the existence of ferromagnetism below Tc is also 
observed. In III-V DMS samples, the hysteresis curves 
are clearly observed as well. However, even at rather 
large fields H, M{H) does not saturate, but continues 
to increase with increasing magnetic field. This featucej 
has been attributed to a "paramagnetic" component ,E3 
and it obviously corresponds to the weakly-interacting 
component of nearly-free spins of our simplified two- 
component model. In fact, we have generated M{H,T) 
curves within the mean-field approximation, and these 
features are clearly present (see Fig. |l^), in qualitative 
agreement with measurements. 

E. Effect of other interactions 

We have investigated in detail the effect of the on- 
site disorder term idi^ and of the on-site Coulomb re- 
pulsion U elsewhere.Eil The on-site disorder is due to the 
Coulomb potential created by the charge impurities re- 
sponsible for compensation (such as AS++ antisites). We 
have considered two extreme possibilities. In the first 
case, we assume that these potentials are completely un- 
correlated, and model them by choosing random values 
for u{i) within an interval [— The estimate W w 
1 Ry is obtained foUojwing standard considerations for 
doped semiconductors In the second case, we attempt 
a simple modeling of the effect of As antisites. We choose 
random positions for these As defects on the Ga sublat- 
tice and identify their two nearest neighbor Mn sites. 
Each such As impurity has an effective charge +2e, and 
therefore will contribute an on-site Coulomb potential 
+2e'^/er at a Mn impurity site which is at a distance 
r from it. However, since the Mn ions also have effec- 
tive ionic charge — e, the As potential is screened (par- 
tially compensated) by the potential of the Mn impurities 
nearby it. Therefore, we assume that each As antisite 
only contributes to the on-site potential u{i) of its two 
nearest Mn neighbors, with the contribution to the other 
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Mn sites being screened out by the contribution of these 
two nearest Mn sites. The presence of the charged impu- 
rities responsible for compensation increases the amount 
of disorder (inhomogeneity) in the system, since the holes 
will avoid the regions were these defects are located. 
Thus, one might assume that u{i) ^ will lead to a fur- 
ther increase of Tc- However, in fact we find a decrease 
of Tc for these mpdels of compensation, especially for 
the second model£2l This is a consequence of the fact 
that due to the presence of nearby As antisites, holes 
now avoid some Mn sites that would otherwise be part 
of dense clusters. Thus, the system effectively moves to- 
wards the more homogeneous regime, with lower Tq- An 
opposite effect is provided by the on-site electron-electron 
Coulomb repulsion, the presence of which was found to 
lead to an increase of Tc, since it aids in the splitting of 
the up and down spin bapds, favoring spin polarization 
at higher temperatures.Eil 

A quantitative determination of the effects of these 
types of interactions will have to wait until more details 
are known about the compensation processes. A theory 
that properly and self-consistently describes the screen- 
ing processes is also necessary. 

V. CONCLUDING REMARKS 

In this paper, we have discussed the behavior of a 
model of DMS in the low density regime, based on a 
simple tight-binding hydrogenic model of the impurity 
band. Such a model takes into account, at the very 
outset, the inherent disorder present in the experimen- 
tal system, namely the random position of the dopants. 
Other modelg3 start from an electron gas exhibiting the 
translational symmetry of the host lattice and ignore the 
disorder of the alloy system. While the latter may be the 
appropriate starting point for the high carrier density 
regime, it does not allow for a metal-insulator transition, 
and consequently misses the unusual transport and mag- 
netic behavior associated even with metallic systems in 
the vicinity of such a transition. In contrast, our model 
starts from the low density insulating side, and at least 
for conventional doped semiconductors, has been found 
to be applicable to densities ub to a factor of 3 above 
the metal-insulator transitionEO 

For the case of II- VI DMS, we have restricted our- 
selves to low densities corresponding to the insulating 
phase, for the case of a half-filled band i.e., no com- 
pensation. In this limit, a spin-only description of the 
bound carriers is appropriate. [We note, however, that 
such a spin only description has been very successful for 
the low temperature thermodynamic and magnetic prop- 
erties in conventional doped semiconductors, both un- 
compensated and compensated, for densities up to the 
metal-insulatcuL-jxansitionO and even somewhat into the 
metallic phasecjo provided an itinerant Fermi-liquid like 
second component is added to the description of these 



highly disordered systems]. For III-V DMS, where large 
compensation is found to be experimentally present, pre- 
sumably due to antisite defects, we have adopted a full 
fermionic description of the carriers. Such a model al- 
lows for both an insulating and a metallic phase. How- 
ever, as explained in the body of the paper, the model 
we have studied is simplified, and neglects several terms 
in the full many-body Hamiltonian describing these com- 
plicated materials. 

Despite the rather different model descriptions (spin 
vs. fermion) for the two cases, as well as methods of 
solution (Monte Carlo vs. Mean Field Approach), we 
find a remarkable similarity in the qualitative predictions 
concerning the magnetic and thermodynamic proper- 
ties. Most striking are the unusual magnetization curves 
Af(T), with linear to concave upwards shape over much 
of the ferromagnetic region, in striking contrast to con- 
ventional uniform ferromagnets. This appears to be a 
combined result of low carrier density and strong disor- 
der. As a consequence, the ferromagnetic transition has 
percolation like characteristics, with only a small fraction 
of the material carrying the bulk of the ferromagnetism 
around Tc. The remaining portion of the material orders 
gradually as the temperature is lowered, and unlike in 
most conventional ferromagnets, saturation magnetiza- 
tion is not reached until well below Tc. Such an inhomo- 
gencous magnetization results in unusual susceptibility 
and specific heat in the low temperature ordered phase, 
and would imply substantial inhomogeneities in the local 
field at Mn sites, which could be probed, e.g., by NMR 
measurements. Unusual hysteresis curves in M{H) be- 
low Tc are also implied, with saturation occurring well 
beyond where the loops close. We have checked for the 
case of the III-V DMS that within a simple impurity 
band description, these effects are robust.EJ However, as 
the carrier density is increased (by reducing the compen- 
sation, or raising the Mn concentration), the anomalous 
shape of Af(T) becomes less prominent: M(T) assumes 
the convex upward shape of conventional uniform ferro- 
magnets, and the ensuing unusual properties discussed 
above gradually fade away. 

In contrast to the qualitative shape of the magnetiza- 
tion curves and the ensuing thermodynamic and mag- 
netic behavior, the actual transition temperatures of the 
two systems are known to be rather different ( from 
a few de^ees. JtelvinB for the II- VI DMS, to several 
hundredsiNlsa-Ell for the III-V DMS ). Certainly one rea- 
son for this difference is the increased weight of the hole 
wave- function at the cation (II/III) site where the Mn 
spin resides in the III-V semiconductors relative to the II- 
VI semiconductors, as may be_scen from a tight-binding 
description of valence bandsH of zinc-blende structure 
semiconductors. However, an additional reason within 
an impurity band description of carriers, is that in III- 
V DMS the Mn sites are centers of the impurity wave- 
functions, while for II- VI DMS, the carrier impurity sites 
are distinct from the Mn. Consequently, the Mn sites 
see a lower amplitude of the carrier wave-function, and 
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thus a lower effective exchange coupling in the II- VI. This 
peaking of the impurity wave-function at the Mn site in 
the III-V based DMS, leads to a further enhancement of 
their Tc vis-a-vis the II- VI based DMS. 

For both the II- VI and the III-V DMS, we find that 
Tc is enhanced by disorder. This can be understood by 
recognizing that in a heavily disordered system, nature is 
able to create global ordering by finding the tortuous per- 
colative pathway necessary when the average coordin^ 
tion number is much below that of any uniform lattice.Ej 
In III-V, the large compensation adds an additional de- 
gree of freedom to the carriers, in the choice of amplitudes 
on different sites, which again leads to wide variation in 
the effective fields at different sites, and implies a percola- 
tive aspect to the magnetic ordering transition. In mean 
field, we find the enhancement of Tc to be quite large; 
however, preliminary Monte Carlo results suggest lower 
effects of jdisorder on Tc than given by the mean-field 
approach. EHl 

Disorder effects on the electronic wave-functions 
will lead to significant transport anomalies, especially 
near the metaj-, insulator transition, as has been seen 
experimentally.El It will also likely affect the nature and 
amount of magnetic scattering of carriers injected into 
the system. While the naive expectation is that disor- 
der should increase spin-flip scattering, it may be signif- 
icantly reduced for carriers near the Fermi level. This is 
because we find that these states have large amplitudes 
along the percolating backbone of the system, where 
the Mn moments are magnetized well above the aver- 
age global magnetization. In this regime, many standard 
models devised for translationally invariant systems {e.g., 
relationship of anomalous Hall effect to bulk magneti- 
zation) may not be applicable, and such interpretations 
should be used with care. 

One other approximation inherent in our work is the 
assumption that electron and hole doping give rise to 
Hamiltonians that are qualitatively similar, though quan- 
titatively different (holes have angular momentum 3/2, 
while electxons have spin 1/2). This is what is found for 
free holesS although the more complicated anisotropic 
wave- functions for holes leads to quantitative differences, 
qualitatively the results are similar, in that both systems 
lead to ferromagnetic ordering. Recently, however, sinpa 
the Spintronics 2001 conference, it has been proposedEil 
that spin-orbit coupling can lead to effective spin-spin 
couplings that are anisotropic in spin space, and the po- 
sitional disorder effectively leads to random anisotropy. 
This could lead to frustration effects not present in our 
model, and if true, would need to be put in for hole doped 
systems to achieve full understanding of magnetic order- 
ing and carrier transport in DMS systems. 

Finally, we discuss the applicability of our results to 
actual III-V DMS in the regime of large Tc- While our 
model is based on the insulating, low density limit, how 
many of its features persist into the metallic phase at 
higher densities and temperatures, is dependent on the 
nature of the filled electronic states at temperatures Tc 



and below. In the model we have studied, the host va- 
lence bend is completely neglected, and its inclusion is 
not expected to lead to qualitative changes, because it 
lies several hundred meV above the Fermi level. How- 
ever, this is a consequence of an impurity band with a 
density of states that is characteristic of a bandwidth of 
order hundred meV also. Should the impurity band be- 
come much broader in the actual system due to effects 
we have left out, it will likely merge into the host valence 
band, and the states will be strongly mixed. Neverthe- 
less, the occupied states for small filling (low Mn density 
and large compensation) would have significant effects 
of disorder. This, in turn, implies that the anomalous 
behavior exhibited by our model would be present, but 
with lower magnitude than shown by our calculations. 
The clearest signature of these would likely come from 
local probes, which would be able to determine the dis- 
tribution of local fields, and hence local density of states 
at various sites. Such input into phenomenological mod- 
els should provide a fruitful avenue for a more in-depth 
study of the fascinating world of real Diluted Magnetic 
Semiconductors, which offer both a significant promise in 
terms of their applications in spintronics, and a challenge 
in terms of their fundamental understanding. 
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FIG. 1. Dependence of reduced magnetization M/Mo 
upon reduced temperature T/Tc. Curves are slightly differ- 
ent for different values of the quantum spin S, however they 
all have a convex upward shape. The solid circles represent 
typical experimental data, for Gd [S ~ |), Fe [S ~ 1) and Ni 
(S ^ i). (From Stanleyll3). 
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FIG. 2. Schematic comparison of typical experimental 
measurements of temperature dependence of magnetization, 
specific heat and susceptibility for a Heisenberg ferromagnet 
(such as EuS) with the predictions of the Weiss mean-field 
theory. Note that the curve for inversepsusceptibility 1/xt is 
shown only for T > Tc. (From Stanley^). 
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FIG. 3. Schematic drawing of the percolation hmit for 
a disordered collection of BMPs. As the temperature is low- 
ered, the size of each BMP increases and a percolated network 
appears below Tc- Just below Tc, only a small fraction of the 
spins belong to the percolated network and sustain the bulk 
magnetization of the sample. The large majority of spins 
is outside the percolated network and behave like quasi-free 
(non-interacting) spins. 
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FIG. 4. Magnetization per Mn spin as a function of 
temperature, in a II- VI DMS, for classical/discrete spin 
model (left/right panel). Results are shown for samples with 
A'' = 256, 864 and 2048 Mn spins. Finite size-scaling analysis 
finds a critical temperature Tc = 0.014Jo.tJ The magneti- 
zation curves are very unlike the conventional ferromagnet 
magnetization curve shown in Fig. |l| 
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FIG. 5. Specific heat per Mn spin as a function of temper- 
ature, for classical (empty squares) and discrete (full squares) 
spin models. Systems with N = 2048 Mn spins were used in 
both cases. While the discrete model recaptures the proper 
limit Cv — > as r — > 0, the peak in the specific heat is 
well below the critical temperature {Tc = 0.014Jo for these 
parameters), unlike in conventional ferromagnets, where the 
peak in the specific heat is at Tc (see Fig. H). 
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FIG. 6. Magnetization per Mn spin as a function of teip, 
perature, in a III-V DMS, using Monte Carlo simulations.E3 
Curves correspond to j: = 0.01, and relative hole to Mn con- 
centrations p = 10 and 30%. 
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FIG. 7. Average magnetization of the Mn spin {Sm„ > 0) 
and of the charge carrier spins Sh < 0) as a function of temper- 
ature, in a III-V DMS. Curves correspond to j: = 0.0093 and 
p = 10 (full line) and 30% (dotted iliae) , and were obtained 
using the mean-field approximationEil 




FIG. 9. Histogram of effective couplings pi = Jeff{i)/J of 
different Mn spins at ksT/J = 0.01, for x=0.01 and relative 
hole to Mn concentration p — 10%. |— This distribution was 
found using Monte Carlo simulations E3 




FIG. 8. The average Mn spin Smti and average spin per 
hole Sh for doping concentration x — 0.00926 and p = 10%. 
In increasing order of Tc, the curves correspond to ordered, 
weakly disordered, moderately disordered and completely 
random distributions of Mn (see text). 




FIG. 10. Hysteresis curves obtained within the mean-field 
approximation for one disorder realization corresponding to 
a Mn concentration x — 0.03 and hole to Mn ratio p — 10% 
(corresponding to a critical temperature ksTc / J = 0.85). 
Averages over several disorder realizations are needed to ob- 
tain smooth curves. 
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